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Abstract 

We describe a procedure for determination of freeze-out parameters of heavy-ion collisions 
through direct comparisons between experimentally measured higher order cumulants of charge 
fluctuations and first principle (lattice) QCD calculations. 



1. Introduction 

The QCD critical point is a conjectured second order phase transition point in the temperature 
(T) and baryon chemical potential (/ig) plane of the QCD phase diagram. Being a second order 
phase transition point it will be associated with large correlation lengths, which in turn will be 
manifested in characteristic large fluctuations. In this vein, one of the major focus in search 
for the QCD critical point in the Relativistic Heavy-Ion Collider (RHIC)'s Beam Energy Scan 
(BES) program has been concentrated on measuring higher order cumulants of conserved charge 
fluctuations, such as net baryon (B), net electric charge (Q) and net strangeness (S). 

Experimentally measured hadronic observables in Heavy-Ion Collisions (HIC) characterize 
the freeze-out condition. The success of statistical hadronization models, based on the ther- 
mal hadron spectrum of the Hadrom Resonance Gas (HRG) model, in fitting the experimentally 
measured hadron yields suggests that freeze-out conditions in HIC can be described by equi- 
librium thermodynamics characterized by freeze-out temperatures (jf ) and chemical potentials 
(p. B ,/j.g,f/ s ). On the other hand, the QCD critical point will be located along the QCD transi- 
tion/crossover line in the T - jjg plane. Thus in order to observe signatures of the critical point 
in RHIC BES, the freeze-out of the conserved charge fluctuations must happen at some (T^,/i B ) 
close to the QCD transition line in the T - fig plane. While for moderate values of fig the QCD 
transition line in the T -fi B plane is known from first principle Lattice QCD (LQCD) calculations 
|Q][2), till now there is no first principle QCD determination of freeze-out parameters associated 
with the observables related to the conserved charge fluctuations. 

Here we outline a procedure for the determination of the freeze-out parameters (7^ , /j. B , /j.q, ^ s ) 
via direct comparison between first principle LQCD calculations and the experimentally mea- 
sured cumulants of charge fluctuations. While such a procedure will tell us whether the exper- 
imentally measured fluctuations can indeed be described by equilibrium thermodynamics, the 
freeze-out parameters obtained using this procedure do not necessarily correspond to the chem- 
ical freeze-out parameters. Chemical freeze-out indicate irrelevance of inelastic scatterings, but 
the freeze-out of the conserve charge fluctuations takes place through diffusion processes. It is 
likely that hadronic final state interactions are important for the freeze-out of charge fluctuations 
and it is a-priori not clear that the charge fluctuations freeze-out at the chemical freeze-out point. 
Thus, in order to make ab-initio parameter-free theoretical predictions relevant for the RHIC BES 
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Figure 1: (a) Top: Continuum extrapolated LO in fig contribution for the electric charge chemical potential as a function 
of temperature. Bottom: NLO contribution, in units of the LO contribution, as a function of temperature for the electric 
charge chemical potential, (b) Same as the previous panel, but for the strangeness chemical potential, (c) Electric charge 
(top) and strangeness (bottom) chemical potential as a function of fig for the relevant temperature range T = 150 - 170 
MeV. 



it is crucial to know the freeze-out parameters associated with the conserved charge fluctuations. 
Once these freeze-out parameters are known from the comparison of the lower order cumulants 
only then higher order cumulants can be predicted from LQCD calculations along this freeze-out 
line in a completely parameter-free manner. Subsequent comparisons with the experimental data 
will clarify to what extent the higher order cumulants contain non-critical/critical signatures. 

The LQCD computations that will be presented here were performed using 2+1 flavor Highly 
Improved Staggered Quarks (HISQ) with a physical strange quark mass and light quark masses 
corresponding to the Goldstone pion mass of 160 MeV. Three different lattice spacings corre- 
sponding to the temporal extents N T = 6,8, 12 were used for the calculations. Further details 
concerning the LQCD calculations as can be found in Refs. [UEl. More detailed discussions 
regarding the determination of the freeze-out parameters can be found in Ref. J4J. 



2. Electric charge and strangeness chemical potentials 

By using the constraints of initial strangeness neutrality and initial iso-spin asymmetry of 
the colliding nuclei of HIC the electric charge and strangeness chemical potentials, fig and /j. s 
respectively, can be treated as dependent parameters of T and hb- Assuming spatial homogeneity, 
the initial strangeness neutrality gives (ns ) = and the initial iso-spin asymmetry of the colliding 
nuclei leads to (hq) = r(ng). Here, «x denotes the density of the corresponding net conserved 
charge X and r - N P /(N P + N n ) is the ratio of the total number of protons to the total number of 
protons and neutrons of the initially colliding nuclei. We choose to work with r = 0.4, a good 
approximation for Au - Au as well as Pb - Pb collisions. By making Taylor expansions of (rix) 
in (jiB,fiQ,fJ-s) U P t° 0(i/ x ) and imposing the above two constraints one can write down the {1q 
and fj.s in terms of the other two independent parameters: hq(T,hb) = q\(J)[iB + Q3(T)lu B and 



In Fig. 1(a) we show LQCD results for the Leading Order (LO) contribution q\(T) (top) and 



the Next-to-Leading Order (NLO) contribution qi(T) (bottom) to fiQ. Similar contributions for 



the jig are shown in Fig. 1(b) These results show that the NLO contributions are less than 10% 
and well under control for a wide range of the baryon chemical potential fis < 200 MeV, i.e. for 
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Figure 2: LQCD results for the thermometer (a) and the baryometer R® (b). (c) LQCD results for the double ratio 
Rp/Rp as a function of Hb/T ■ 



19.6 GeV. In Fig. 1(c) we show hq{T,[i b ) (top) and fis(T,fiB) 



RHIC energies down to ■sjspft 
(bottom) as a function of fig for the relevant temperature range T = 150-170 MeV. Comparisons 
with the HRG model suggest that the LQCD results differ by 10-15%. 



3. Freeze-out temperature and baryon chemical potential 

As fig and fj. s are known as a function of T and fi B all the cumulants of conserved charge 
fluctuations can be now be expressed only in terms of T and /i^. Although in RHIC both net 
proton and net electric charge fluctuations are being measured J5] [6] [7), in LQCD only fluc- 
tuations of conserved quantities are accessible. Since the net proton fluctuations may not be 
quantitatively equal to the net baryon number fluctuations, it is safer to work with the net electric 
charge fluctuations. We propose to look at two different volume independent ratios formed out 
of the three lowest order cumulants, mean (Mq), variance (<x G ) and skewness (S q), of the charge 
fluctuations — 
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While the cumulants of the net charge distribution Mq, ctq, S q are being measured in RHIC BES 
B|6] EL the generalized higher order charge susceptibilities Xn(T,/J-B) — d"lnZ,/d(jj.B/T)" can 
be calculated using LQCD via Taylor expansions Xn(T>MB) = 'Lk=oX^+^T' Q)Q*b/T)" /k\. In a 
similar manner we can Taylor expand the ratios R® { and R® 2 themselves up to NLO in fig as 



shown in Eq. (|la|i and Eq. (|lb|) respectively. Note that, in the LO the ratio R^ is independent 



of /j B but the LO term for the ratio R® 2 is proportional to Thus, R® can be used as the 
thermometer to determine and R® 2 can used as the baryometer to fix [J B . 

In Fig. 2(a) we show the LQCD results for the ratio R^ . The estimated NLO correction to this 
ratio is at most 1 0% for / T ~ 1 over the whole T range. The broader band in the figure depicts 
the range of R? x for fis/T = 1 including the NLO correction and the thinner band depicts the LO 
results, i.e. for pt B - 0. The LQCD calculations for the thermometer R^ shows a characteristic 



T dependence and large deviations from the HRG results within the relevant temperature range 
T = 150 - 170 MeV. For example, by comparing with these QCD results an experimentally 
measured value of > 2 will indicate a freeze-out temperature T-f < 150 MeV, ~ 1.5 will 
give T f * 160 MeV and < 1 will mean T f > 170 MeV. 

After fixing from the thermometer Rr x we can use the baryometer R® 2 to determine the 
freeze-out baryon chemical potential fJ B . Fig. 2(b) shows LQCD results for the ratio R® 2 as 
a function of fig in the relevant temperatures interval T = 150 - 170 MeV using continuum 
extrapolated LO results and adding contributions up to NLO in fi B . In this temperature range the 
NLO corrections are well under control, less than 10%, for fi B < 200 MeV, i.e. for RHIC BES 
energies of ysMV ^ 19.6 GeV. By comparing the experimentally measured values of R® 2 with 
these first principle QCD calculations one can determine fi B . As a concrete example, choosing 
Tf = 160 MeV, an experimental value of R® 2 = 0.01 - 0.02 will give fJ B IT s = 0.1 - 0.2, R.f 2 = 
0.03-0.04 will suggest f/jT^ = 0.3-0.4 andfif 2 = 0.05-0.08 will indicate f/ B /Tf = 0.5-0.8. 



4. Thermodynamic consistency 

Such a determination of the freeze-out parameters for a given beam energy will help us in 
understanding whether at that beam energy higher order cumulants of conserved charge fluctu- 
ations can be consistently described within the framework of equilibrium thermodynamics. If 
the experimentally measured fluctuations are indeed described by equilibrium thermodynamics 
characterized by unique values of temperature and chemical potential then other volume inde- 
pendent ratios of cumulants must also have uniqu e values as predicted by (L)QCD calculations 



at the same {T? ,fi B ). As an example, in Fig. 2(c) we show the LQCD predictions for the double 
ratio Rf 2 /^i2 as a function of fisIT . For a pre-determined values of the freeze-out parameters 
and ft B this ratio has a unique value consistent with equilibrium thermodynamics. If the ex- 
perimentally measured observables contain the correct physics of equilibrium thermodynamic 
fluctuations then they must agree with this (L)QCD prediction. 
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